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Abstract. Within a modified Poisson-Boltzmann theory we study the effect of 
Bjerrum pairs on the typical length scale l/S over which electric fields are screened in 
electrolyte solutions, taking into account a simple association-dissociation equilibrium 
between free ions and Bjerrum pairs. At low densities of Bjerrum pairs, this length scale 
is well approximated by the Debye length 1/k oc with ps the free ion density. 

At high densities of Bjerrum pairs, however, we find 1/k cx which is significantly 
larger than 1/k due to the enhanced effective permittivity of the electrolyte, caused 
by the polarization of Bjerrum pairs. We argue that this mechanism may explain 
the recently observed anomalously large colloid-free zones between an oil-dispersed 
colloidal crystal and a colloidal monolayer at the oil-water interface. 
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1. Introduction 

Phase separation and criticality in electrolyte solutions, following from Debye-Hiickel 
(DH) theory [1], was extensively studied by Fisher and Levin [2]. The calculated 
critical density and temperature differed by only ~10 % from the results of Monte 
Carlo simulations by Panagiotopoulos [3], showing DH-theory to be a reliable basis to 
describe some basic features of electrolyte solutions. Fisher and Levin extended the 
original DH-theory by inclusion of Bjerrum pairs, following the ideas of Bjerrum [1] 
that plus and minus ions can form neutral pairs. These Bjerrum pairs or dipoles are 
considered as a separate particle component, and reduce, within the restricted primitive 
model, the total number of free ions. The results of the extended DH-theory agreed 
remarkably with simulation results, especially when couplings between the dipoles and 
the ions, as well as hard-core repulsions were taken into account [2]. Here, we apply 
the same ideas of Bjerrum to describe screening effects in low-dielectric solvents ('oils') 
by means of a modified Poisson-Boltzmann theory. In these solutions we expect strong 
correlations between the ions, since the energy gain of bringing two oppositely charged 
ions at contact can exceed the thermal energy considerably. 

We will consider systems of free ions and dipolar particles, similar to the study of 
e.g. [5J, and find that electric fields are screened over a typical length scale 1/k that can 
be significantly larger than the Debye length 1/k (based on the ionic strength), at least 
for large dipole densities. We predict that these densities are to be expected in low- 
dielectric solvents by considering a simple association-dissociation equilibrium between 
the free and bound ions |6]. Recently, strong electrostatic repulsions were observed by 
Leunissen et al. [7] in low-dielectric solvents (4 < e < 10), sometimes extending over 
a length scale beyond 100 /im. The Debye length 1/k^^ was found to be only ~ 4 
fim, calculated from conductivity measurements (the density of free ions). The analysis 
presented in this paper provides a possible explanation and a reason for the quantitative 
differences between these experiments [7] and the theory in [8], where Bjerrum pairs were 
not taken into account. 

In section |2] a reaction equilibrium between free and paired ions will be discussed. 
Parameter space will be divided into regions where bound ions outnumber the free 
ions, and vice versa, as was already calculated, e.g. in [6]. In section |3] the effect of the 
dipole density on the screening length will be analyzed. First we extend Gouy-Chapman 
theory [9] by including dipoles and calculate the effective screening length 1/k, similar 
to the calculations in [5j. Remarkably, we find 1/k oc ^/p^ at high dipole densities, in 
contrast to 1/k oc 1/ s/pl, where ps is the density of free ions. Finally we review and 
extend the theory presented in |8] and predict a larger effective screening length due to 
dipoles. The considered densities of free ions are well below the critical value, i.e. phase 
equilibrium is not considered, even though the temperatures of interest are close to the 
critical temperature. 
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2. Bjerrum pairs 

First we consider a three-dimensional bulk electrolyte of monovalent cations and anions, 
at a total density of 2ptot- The ions may form pairs or remain free; the number density 
of dipolar Bjerrum pairs is pd and the number densities of free ions are p+ = p_ = ps. 
In terms of dimensionless densities t]^ = PxCr^, where a is the common diameter of the 
ions, the total density ?7tot of ions (of one type) is 

Vtot =r]d + Vs- (1) 

The strength of the electrostatic interactions in the solvent is reflected by the Bjerrum 
length, which is the length at which the bare Coulomb interaction between two 
monovalent ions is exactly 1 ksT, 

- 4^' PI 

where e is the elementary charge, and e the dielectric constant of the medium. For water 
at room temperature, this length is only 0.71 nm, for apolar solvents it measures up 
to several tens of nanometers, and in vacuum it is ~ 57 nm. The Coulomb interaction 
between two ions can hence be written in terms of the Bjerrum length, 

w - ^- - 4- 

KbT r I 

where the +, — refer to equal charged particles, and oppositely charged particles, 
respectively, r is the distance between the particles, and I is a dimensionless distance. 
The dimensionless equilibrium constant K of the reaction of free ions that bind into 
paired ions is defined by 

A- = !?±!h = !l = £!e.p(^). (4) 

Vd Vd \kBTJ 

where A is the ionic Debroglie wavelength, and where AG is the free energy of a bound 
pair of ions (being separated by a distance a < r < Xb), with an associated Coulombic 
binding energy Vq = —k-BT/l. It will be convenient to introduce the dimensionless 
temperature T* = (t/Ab, such that can be expressed in terms of an internal partition 
function 

1 

/ 1 \ 3 /■ /I 



K-'=Arr[-) J dl exp(^-j, (5) 

T* 

as already postulated by Bjerrum It can easily be checked that T* ~ 1 for typical 
ions such as Na"*" and Cl~ in water at room temperature, and T* < 0.2 in oils with e < 15. 
Figure [T] shows the relation between the equilibrium constant K and the dimensionless 
temperature T*. With increasing temperature K rises, thereby lowering the tendency 
to form pairs according to (jl]). Using this can be further quantified by relating ris to 
Vtot as 

2 

Vtot =Vs + ^, (6) 
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which yields 1]^ and t]^ as a function of r^tot and K{T*) as 
??tot 2?7tot W K / ?7tot' 



* 

T 
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Figure 1: The dimensionless equiUbrium constant Ka^ = K as a function of the dimensionless 
temperature T* = ct/Xb- The hne also marks the points where the density of paired ions equals 
the density of free cations/anions for given concentration of free ions (see right vertical axis), 
for an ionic diameter a = 0.5 nm. 

From (jlj) it can be easily seen that ii K = 7]^, then ?7s = 77^. If the left vertical axis 
of figure [1] is read as the dimensionless density rjs of free ions, the curve K{T*) thus 
separates the parameter regime where dipoles dominate (?7d/''?s > 1)? from the regime 
where free ions dominate {fjd/Vs < !)• The right vertical axis converts the corresponding 
?7s = K{T*) to the molar density ps for the typical choice a = 0.5 nm. We have already 
seen that T* ^ 1 for typical ions such as Na^ and Cl~ in water at room temperature, 
and that T* < 0.2 in oils with e < 15. Figure [1] therefore illustrates, for instance, that 
an electrolyte with a millimolar ionic strength, ps — 1 mM, is dominated by free ions 
for T* > 0.15 and by dipoles for T* < 0.15. For nanomolar concentrations, Ps = 1 nM, 
the crossover is at T* ~ 0.05. For later reference we also consider the mean separation 
a between the ions in a pair. Assigning a statistical weight oc exp(AB/r) to a pair at 
separation r, one finds 

1 

= = a-KA.[-^y j d/ exp (i), (8) 

which can straightforwardly be evaluated numerically. The result is shown in figure [21 
as a function of T*, for several upper bounds of the integration domain to check the 
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dependence of a and K on the definition of a Bjerrum pair, being two oppositely charged 
particles separated by a distance < Ab- At low temperatures T* < 0.05 the values do 
not depend on the precise definition, because the deep potential well of > 20 k^T at 
contact dominates th( 

2.5 



2 
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lb 
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Figure 2: The mean ion separation of the Bjerrum pairs a (in units of the ionic diameter 
a) as a function of the dimensionless temperature T* = a/XB, for several definitions of the 
maximum ion separation that is still called a Bjerrum pair. At low temperatures T* < 0.05 
the separation does not depend on the precise definition, because the deep potential well of 
> 20 k^T at contact dominates the partition sum of the pair. 




3. Effective screening length 

We now consider a system of monovalent ions near a charged plate at 2; = 0, where the 
2;-axis is perpendicular to the plate, and the ion density is ps at 2; 00. We expect 
that the charge of the plate is screened by an oppositely charged ionic cloud of net 
charge, generated by the ions. By the Poisson-Boltzmann equation, solved by Gouy 
and Chapman [S] , the typical width of this cloud (the double layer) can be found. This 
length scale is also known as the Debye length 1/k = 1 /v^SttAbPs- We now extend 
Gouy-Chapman theory |9j by the inclusion of an additional particle species of dipoles 
with number density pd in the bulk far from the plate. Following the derivation of |5j, 
we find the Poisson-Boltzmann equation for z > 

<P"iz) = K\mh<Piz) - n^f^a^^[g(a<P'izm (9) 

where a prime denotes a derivative with respect to z, and / k-QT(f){z) /e is the electrostatic 
potential. The function Q{u) = cosh(u)/u — sinh(u)/M^ can accurately be approximated 
by a first order expansion since the mean separation cr = (9(1) 

nm, and the electric field (f)'{z) = 0{1) pm^^ for the systems of our interest. The 
PB-equation then reduces to 

(p"{z) = R'^smh^{z), (10) 
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where 

= T , (11) 

with a = p^a^/Gps. The presence of dipoles thus increases the screening length 
significantly as soon as aK^ = 0{1), which can only be obtained at high ionic strength 
in low dielectric media, such that pd/Ps is large. Equivalently, one can also consider 
the dielectric constant to be effectively changed by the presence of the dipoles. Writing 
R = 8ne'^ps/{ekBT), with e the effective dielectric constant gives with (fTTj) that 

. = . + i^^!^. (12) 

3 

where is the Bjerrum length in vacuum. The molar density of pairs has to be 
large enough for a significant change in the effective dielectric constant. For typical ion 
diameters of a few Angstrom one needs pd ^ 10 mM for e > 2e. 

In (l9l)-( fT2l) we treated ps and pd as independent densities, whereas in (jll) we related 
them directly through an equilibrium reaction. Using (jl]) we find a = risa'^/{QK) such 
that the limit an"^ S> 1 gives 

= y/a oc ^/p^, (13) 

which is a remarkable dependency, since in the absence of dipoles the Debye length 
scales as k,~^ oc 1/ ^fpl. The full dependence of on follows from (fTTj) and (jlj) as 

{^2 
-=1 ,'«"'«^; (14) 

with = T*a^a'^ /(ASnK). Relation (1141) is plotted, in a conveniently scaled fashion, 
in figure [3^, revealing a minimum at = 2^ A that separates the dilute limit (where 
K ~ k) from the dense limit (where A^k). In figure [3b the ratio k ^/k ^, which 

equals a/ {AkY + 1 from f|T4|) . is plotted as a function of T* for pd = 0.1, 1, 10 M, setting 
a = a = 0.5 nm. Figure [3] reveals a strong modification of the effective screening length 
(by factors of 2 - 10) provided pd > 1 M and T* < 0.1 (or e < 10). The key question is 
therefore if such conditions are experimentally attainable; the answer will be provided 
in the next section. 



4. Physically relevant regime 

From and figure [3] we conclude that the effective screening length is significantly 
larger than the Debye length, if the density of dipoles is of the order of 1 M. We do not 
expect that these densities can be reached in polar solvents (water), where T* is high, 
and hence K is high, such that free ions dominate according to ([7]) and figure [H The 
equilibrium constant K can be very low in low-polar solvents (oil), such that there are 
many more Bjerrum-pairs than free ions. However, free ions also dissolve worse in oil 
than in water, and the question is whether or not enough ionic strength is possible to 
produce the Bjerrum pairs at all. In order to get an estimation of the parameter regime 
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Figure 3: The effective screening lengtfi as a function of tfie Debye length k ^ 

^g^j^ I J: ;„„! J. J. ^^* „„ „ „f < piT[ \ 
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(a) The (dimensionless) effective 
screening length as a function of 
the Debye length showing the 
asymptotic regimes — 
and ~ kA^, separated by a 

minimum at k^^ = 2* A (see text). 
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(b) The effective screening length (in units 
of the Debye length) as a function of 
dimensionless temperature T* for several 
dipole densities — 0.1, 1, 10 M, at a mean 
separation and particle diameter a ^ a = 0.5 
nm. The upper axis shows the conversion from 
r* to e at cr = 0.5 nm. The effective screening 
length can be up to a factor 0(10) higher than 
the Debye length at high dipole densities. 



where inflation of the screening length could take place, we calculate the minimal free 
ion density p™™cr'^ = r/™™ at which k/k > 1 to a significant degree. From (jlj) it can be 
found that a/t^ > 1 implies 



where the dimensionless f]s,T*, and K were defined in the previous section. This 
condition provides a lower bound for the ion concentrations in oil where a significant 
effect from pairs can be expected. If we consider the oil to be in contact with a water 
reservoir, and assume the ions to partition between the two phases due to a difference 
in Born self-energy p!2| [T3| [HI [T7j , we can estimate the maximum ionic strength p^^^ 
in the oil. Given the maximum ionic strength in water, being about 10 M, we find 

^max _ . 1^ T* j ■ ^ significant increase of k~ due to Bjerrum pairs 

thus requires the existence of a regime where p™™ < ps < p^^^. From a numerical 
analysis, however, we find that p™^"" < p™™ for all T* using reasonable values for a. In 
other words, the required high density of dipoles cannot be reached according to the 
equilibrium constant K defined in ([5]). 

In order to make some further progress, we treat solvation effects in a slightly less 
naive fashion. So far we considered the ions to have an effective radius a± of a few A, 
connected to a Born self energy of several tens of k-^T in oil and less than 1 k^T in 
water. The bare radius of ~ 1 A for small ions such as Na"*" or Cl~ would overestimate 
the self energies and underestimate the solvation of the ions in low polar solvents as 
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found for example in experiments [7|. The actual effective radius is thus larger due to 
hydration shells of water molecules that form a cage around the ion [15]. By assuming 
now that two ions can approach each other up to the bare diameter ^cr of the ions, where 
< ^ < 1, in stead of a, the effective diameter of the ions, the equilibrium constant K is 
found to be much lower. This is visualized by the breaking and forming of the structure 
of water molecules around the ions. Within this speculative picture the equilibrium 
constant of ([5]) is redefined by 

1 

^f^ = 4vr(^)Y dU' exp(i), (16) 

where an explicit energetic and entropic cost of restructuring the layer of surrounding 
water molecules is ignored. A small ^ can lead to much higher densities of dipoles, by 
orders of magnitude compared to ^ = 1. Intuitively one could expect higher order 
clusters to form at such densities. We expect however that higher order clusters 
are increasingly unfavourable. The energy gain by electrostatic arguments has to 
compensate for both the loss of entropy, and the energy needed to restructure the 
surrounding water molecules [T6] . 



Figure 4: The equilibrium constant related to the dimensionless temperature T* = a jXs 
(a), and the densities at which (xk} = 1 (b), for several values of ^. The right vertical axis 
denotes the molar density for the choice cr = 0.5 nm. 
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(a) The equilibrium constant related to 
the dimensionless temperature T* — a/Xs for 
several values of ^. When one focusses on 
the left vertical axis, the lines also mark the 
points where the density of paired ions equals 
the density of free cations/anions for given 
concentration of free ions (see right vertical 
axis), and given an ionic diameter of tr = 0.5 
nm. For decreasing ^ the parameter regime 
where the dipoles dominate increases. 




(b) The ion densities at which aK? = 1, i.e. 
where the effective screening length R^^ = 
2k^^ . For small i^, i.e. a small bare radius 
compared to the effective radius, the screening 
length is increased by the presence of the 
dipoles, already at low salt concentrations. 
The red dashed line denotes the ion density at 
which the salt concentration in water is 10 M, 
i.e. higher salt concentrations are not physical 
in oil. 
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In figure we plot K as a function of T* for several ^ ; for ^ = 1 the curve is 
identical to the curve of figure [H As in figure [1], also the present curves separate the 
high-T* regime dominated by free ions from a low-T* regime dominated by dipoles. 
Reducing the contact distance from a to ^cr is immediately seen to reduce the free ion 
regime: dipoles form already at higher T*. The right vertical axis denotes, again in 
analogy with figure [T|, the free ion concentration at the crossover ps = (using a = 0.5 
nm to convert to molar concentrations). Figure UK thus reveals a lowering of the free 
ion concentration at which dipole formation sets in by orders of magnitude when ^ is 
reduced from 1 to 0.2. Figure shows the (dimensionless and molar) maximum and 
minimum free ion concentrations p™**^ and p™™(0) respectively, as a function of T* and 
several ^. Consistent with our earlier observations we see that p™™ > p™^^ for ^ > 0.6 
in the T*-regime of interest. Interestingly, for ^ < 0.6 a physically attainable regime 
of p™"^ < Ps < p^^^ opens up in which significant dipolar effects are to be expected to 
increase the effective screening length beyond the bare one. For the choice of a = 0.5 
nm used in figure Hb, this implies that Bjerrum pairs play an important role provided 
the bare ion diameter < 0.3 nm. This seems physically reasonable. 

5. Possible observations 

We will now consider the system of [7], consisting of micrometer-sized, strongly 
hydrophobic PMMA particles, dispersed in an oily mixture of cyclohexylbromide and cis- 
decalin in contact with water, containing monovalent ions. A densely packed monolayer 
of colloidal particles was observed at the oil-water interface, and a dilute bulk crystal 
separated by a large colloid- free zone of ~ 100 pm between the bulk crystal and the 
interface. The system was theoretically described in [8] with a model that will be extend 
here by the introduction Bjerrum pairs. We consider strongly hydrophobic colloidal 
particles in oil near a planar oil-water interface, in the presence of monovalent ions. We 
focus on the distribution of particles in the direction perpendicular to the interface. By 
employing the framework of density functional theory we write the grand-potential as a 
functional Q[p, p+, P-, Pd] of the variational density profiles of the colloidal particles p(r), 
the cations p+(r), the anions P-(r), and the dipoles pd(r, s), with s the vector of the 
dipole orientation. It is identically presented in [8] except for the dipole contributions 
(with a subscript d), and given by 




+ 



+ 





drdspd(r, s) (^kBT{\np^{r, s)A^ - 1) + V+(r) + V_(r) + AG(r)) 
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+ kBT j dr(p(r)v^(r/(r)) + iQ(r)0(r)), (17) 

where ^^(r) = 47ra'^p(r)/3 is the colloidal packing fraction, and where the first and 
second line are the ideal-gas grand-potential functionals of the ions and the colloidal 
particles in their external fields, respectively, the third line is the ideal gas free energy of 
the dipoles and the binding free energy, the fourth line a grand canonical contribution 
(for fixed chemical potentials), and the last line describes the hard-core and Coulomb 
interactions, [8j . The chemical potentials of the colloidal particles is represented in terms 
of a reference colloid packing fraction r^o to be discussed below. The total local charge 
number density Q{v) = Zp{r) + p+(r) — p-(r) + J ds[pd{r + s|, s) — pd(i" — s|, s)], with 
s the unit vector denoting the direction of the dipole, and a the mean distance between 
the centers of the ions of a pair (previously found to be a ~ cr in low dielectric media, 
see dHD and figure El). For small a the last term in the expression for Q{v) reduces to 



j ds o-Vpd(r,s) 



■ s. 



We obtain the equilibrium distribution of colloidal particles by minimization of the 
functional (ITTj) with respect to ^^(r), which reduces to ri{z) due to the symmetry of 
the system. Minimizations with respect to the densities p±{z) yield the Boltzmann 
distributions for the free ions, also given in [8] (only we use the slightly different notation 
p{oo) = ps here) and dipole density pd{r,s) 

smh{a)\(p'{z)\ 

Pd{z) = pd — -rrn-Ti — ' (^9) 
a\(f)'{z)\ 

where the expression is integrated over s, and where we used the relations between K, 
Pd and Ps given in earlier sections. Combining these with the Poisson equation yields a 
Poisson-Boltzmann equation (Q from which we can find the electrostatic potential (f){z), 
for the boundary conditions 

lime^0'(2;) = \imeo4>'{z) ; lim 4>'{z) = 0, 

ztO 2J.O z^±oo 

where is the relative permittivity of water, and 60 that of oil. The five equations 
of the five unknowns (the particle densities and the electrostatic potential) are solved 
numerically. 

Figure [5] shows the resulting packing fraction profile ri{z), as well as the electrostatic 
potential (j){z) in the inset, for ^ = 1 and ^ = 0.37. The parameters are the colloidal 
radius a = 1 pm, the colloidal charge Z = 450, the relative permittivity of water = 80, 
and that of oil eo = 5.2, and also the external potentials V{z) and V ± (z) are identical 
to those in [8], i.e. V{z) is a Pieranski potential |T8] with contact angle cos^ = 0.987 
and an oil water interfacial tension of 7ow = 9 mN/m, and V±{z) is based on Born 
self-energy differences in oil and water with ionic radii a±. = 0.3 nm (i.e. a = 0.6 nm). 

The curve ri{z) for ,^ = 1 is virtually identical to the one published in ^ (i.e. the 
effect of the Bjerrum pairs can be ignored completely) and reveals a strong monolayer 
adsorption at 2; ~ a, a colloid-free zone for I ^ z/a < 30, and a colloidal crystal [11] 
with a packing fraction ri{oo) = rjiy = 5 ■ 10~^ at z > 30. Experimentally, however, a 
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Figure 5: The packing fraction profile ri{z) of strongly hydrophobic, oil-dispersed colloidal 
spheres (radius a = 1/im, charge Z = 450) in the vicinity of a planar interface at z = 
between water (z < 0, dielectric constant = 80) and oil {z > 0, Eo = 5.2), for a colloidal 
bulk packing fraction rji, = r]{oo) = 5 x 10~^ of weakly wetting colloidal particles (cos 6 = 0.987) 
and screening length in oil /a = 8. The curves show the influence of Bjerrum pairs on the 
depletion zone. For ^ = 1 the effect of Bjerrum pairs can be ignored {k~^ /a = /a = 8). 
For ^ = 0.37 Bjerrum pairs cannot be ignored and the screening length effectively increases 
up to R~^/a = 20, due to an effective increase of the permittivity, resulting in a long range 
repulsion between the monolayer and the bulk crystal over 100 fim. 

colloid-free zone that extends to z ~ 100 fim was observed [7] for these parameters. 
Comparing this with the curve for ^ = 0.37 in figure yields a much better agreement 
with the experimental observation. We speculate, therefore, that Bjerrum pairing is an 
interesting feature for further study in these oily solvents. 

6. Conclusion 

In this paper we considered the effect of Bjerrum pairs on the screening length, and 
concluded that it can be significantly larger in low-polar media than the Debye length 
that is calculated from the free ion concentration (for example obtained by conductivity 
measurements) and the bare solvent dielectric constant. Due to the coupling of free ions 
and dipoles through an association-dissociation equilibrium, we predict the effective 
screening length to scale as k,~^ oc y/p^ at relatively high salt concentrations, in contrast 
to the scaling oc 1/ ^fpl for the Debye length, where ps is the free ion concentration. 
A large concentration of Bjerrum pairs was found to change the dielectric constant of 
the medium effectively. By a naive treatment of solvation effects of the ions, the required 
dipole concentrations seem to be unattainable for physical parameters, such that the 
effect of Bjerrum pairs could be neglected completely. After making the distinction 
between an effective ionic diameter, due to hydration shells that lower the self-energy, 
and a bare ionic diameter, determining the closest distance between two ions, a regime of 
physical parameters was found where inflation of the screening length could be expected. 
Our results provide a possible explanation for the extremely large colloid-free zone that 
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was observed in recent experiments [7|. Clearly, however, more research is needed to 
investigate this effect. 
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